INFNNA-IV-97/40 
DSFNA-IV-97/40 
hep-ph/9708326| 

August 1997 



Can we extract short— distance 
o\ '• information from B(Kl — > 



< 

(N 



G. D'AmbrosicA G. Isidori and J. Portoles 



a INFN, Sezione di Napoli, Dip. di Scienze Fisiche, Univ. di Napoli, 1-80125 Napoli, Italy 
b INFN, Labor atori Nazionali di Frascati, P.O. Box 13, 1-00044 Frascati, Italy 



> 

(N 
CO 
00 

o 

t^- ■ Abstract 

A new analysis of the long-distance two-photon dispersive amplitude of Kl — ► is 
Qh! presented. We introduce a phenomenological parameterization of the Kl —* 7*7* form 
factor, constrained at low energies by Kl — > r y£ + £~(£ = e, /x) data and at high energies by 
perturbative QCD. Using this form factor we provide a reliable estimate of magnitude and 
relative uncertainty of the two-photon dispersive contribution to Kl — > We finally 

discuss the implications of this analysis for the extraction of short-distance information 
rS ! from B(K L — > 



1 Introduction 



Historically the Ki — > p + p~ decay provided a very important tool for understanding 
the flavour structure of electroweak interactions [[[], |2| and nowadays it still represents 
an interesting window on short-distance dynamics. The amplitude of this process can be 
conveniently decomposed into two distinct parts: a long-distance contribution generated 
by the two-photon intermediate state (Fig. la) and a short-distance part that, within 
the Standard Model, is due to W and Z exchange (Fig. lb). The latter turns out to be 
dominated by the top quark and it is known to the next-to-leading order in QCD 0. 
If we were able to disentangle this contribution from the measured Kl — > branch- 
ing ratio we could extract interesting information on the Cabibbo-Kobayashi-Maskawa 
(CKM) matrix element V t d @- Furthermore, a model-independent determination of the 
short-distance amplitude could be useful to put constraints on possible Standard Model 
extensions ||. 

To fully exploit the potential of Kl — > fi + fi~ in probing short-distance dynamics, it is 
necessary to have a reliable control on its long-distance amplitude. However, the dispersive 
contribution generated by the two-photon intermediate state cannot be calculated in a 
model-independent way and it is subject to various uncertainties || [7], [|, || 10|, The 
purpose of this paper is to re-analyse this contribution, using all available information on 
the Kl — > 7*7* transition and trying to evaluate the error due to the model dependent 
assumptions. We will introduce a new low-energy parameterization of the Kl — > 7*7* 
form factor in terms of two parameters a and (3 measurable from Kl — > ^t r l~ (£ = e, fi) 
and Kl — > e + e~p + fi~. Moreover, we will discuss the matching of this approach with the 
behaviour of the form factor in perturbative QCD. Finally, using our estimate of the two- 
photon dispersive contribution, we will derive new bounds on the CKM parameter p JT2 ] 
and on possible new-physics flavour-changing couplings. 

The plan of the paper is as follows. In Section 2 we briefly discuss the general decom- 
position of the Kl — * p + p~ branching ratio and the main formulae for the bounds on 
short-distance parameters. In Section 3 we introduce our low-energy parameterization of 
the Kl —>■ 7*7* form factor, we discuss the determination of a and j3 and the matching 
with the QCD calculation. Finally, in Section 4, we analyse the numerical results. 



2 Decomposition of B(Kl — > fi + fi ) 

The Kl — > branching ratio can be generally decomposed in the following way 

B(K L -> fi + p-) = \MeA\ 2 + \%mA\ 2 , (1) 

where §leA denotes the dispersive contribution and QmA the absorptive one. The former 
can be rewritten as 

9?CvA tyi-CAlong ~\~ tft&Asjiort ; (2) 
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Figure 1: Long-distance (a) and lowest-order short-distance (b) contributions to Kl 
p + p ~. 



whereas the latter can be determined in a model independent way from the Kl —>■ 77 
branching raticQ 



|Qm.A| 



2 a em m n 



In 



5(^l -> 77) 



\ 



4m?, 

1 - 



(3) 



The recent measurement of B(Kl — > /i + p ) [|J is almost saturated by the value of l^m^Ap 
leaving a very small room for the dispersive contribution M 
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Within the Standard Model the NLO short-distance amplitude can be written as 
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where p — p(l — A/2) [|T5|] and p, A are the usual Wolfenstein parameters (T^|. Using this 
result we can write 



p = 1.2- 



c.cp I ± |9fc&4, 
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m t [m t ) 
170 GeV 



-1.55 



W, 



cb\ 



0.040 



(6) 



where the sign inside the modulus is positive if ?R,eA w and $t,eAi on g interfere destructively 
and |3fJe^uj| > |3?e^ D „ 9 |. In principle the above equation could be used to put both a 

a In principle the absorptive amplitude also receives contributions from intermediate states other than 
two-photons, like the two-pions one, but these are completely negligible 113]. 
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lower and an upper bound on p. However |3?e^4 exp | is compatible with zero and, as we 
will show in the following, the same is true for |3?e^4/ on9 |. Thus the upper bound on p is 
useless since it is above unity. On the other hand, independently of the interference sign 
between $leA w and $teAi on g, we can derive a possibly meaningful lower bound on p 

| ■ (7) 

Beyond the Standard Model we can parameterize new-physics contributions as in |5| , 
introducing a flavour-changing Zds coupling at the tree level. Using the Lagrangian 

cZ NP = j^rT.u,Ai^ L z^ (8) 

Z COS u w ^ • 

we obtain |3?e^.Arp| = 3.7 \tfceUds\- Then, assuming 3teA s hort = ^Anp + $leA w , the most 
conservative bound on |3teL^ s | is given by 

\3leU ds \ < 0.27 max {\3teA exp \ + \$teAi ong \ + |3te4J} • (9) 
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3 The K L — > 7*7* form factor and dleAi ong 

The necessary ingredient for the evaluation of $teAi ong is the construction of a suitable 
Kl — > 7*7* amplitude. Assuming CP conservation, gauge and Lorentz invariance implies 



the following general decomposition |16 



A(K L - 7 *( gi , e 1 ) 7 *( 92 , e 2 )) = ie nupa F(q 2 , q 2 2 ) , (10) 

where F(qf, q 2 ) is a symmetric function under the interchange of q\ and q%, and |.F(0, 0)| 



can be determined by the Kl — > 77 width [17 

'Q4ttT(K l ^77) nV2 



|F(0,0)| 



nr 



(3.51 ±0.05) x lO^GeV" 1 . 



Using (|10|) we obtain 
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The structure of the Kl — ► 7*7* form factor has already been discussed and pa- 
rameterized in different ways in the literature |], |^, |], |TU], TT|. However, all the existing 



analyses use model dependent assumptions and thus suffer from uncontrolled theoretical 
uncertainties. In order to be as model independent as possible and to evaluate the size of 
the theoretical errors, we propose the following low-energy parameterization 

fill <&) = ] ~MM = 1 + a f "IT^T + + @T~2 M-2 2T . ( 14 ) 

j\hx,H2) F(Q;0) \ql-m 2 v ql-m 2 v ) '{q 2 - m 2 v ) (g 2 2 - m 2 v ) 

where a and (3 are arbitrary real parameters and my is conventionally chosen to be the 
p mass. The above expression has at least three interesting features: 

1. It is the most general parameterization compatible with the chiral expansion of the 
K L -> 7*7* amplitude up to 0(p 6 ) ||, [U|. 



2. It includes the poles of the lowest vector meson resonances with arbitrary residues. 

3. The parameters a and (3, expected to be 0(1) by naive dimensional chiral power 
counting, are in principle directly accessible by experiments in Kl — > / ~f£ + £~(£ = 
e,/i) and K L — > e + e~p + p~. 

Clearly the expression (|i~4"D cannot be considered correct for arbitrary values of q\ and 
q\. To be more general we should consider a and (3 as g 2 -dependent couplings. However 
for the purposes of this first analysis (that should be improved in the future along these 
suggestions) we believe it is reasonable to treat a and (3 as constants up to q\ ~ q\ ~ 1 
GeV 2 . Moreover, being just a phenomeno logical description, we do not expect the form 
factor QT4T ) to produce a finite result for the Kl — > p + p~ amplitude. Indeed, using (p~3|) 
and (|14T) we obtain 

UeR{m 2 K ) = -3[ln(A/m ) + 2aln(A/m Q ) + p\a(A/mp)] 

= -3[ln(m /3 /m ) + 2«ln(m /J /m (t )] - 3(1 + 2a + /3)\n(A/m ) , (15) 

where 

m = 140 MeV , m a = 452 MeV , m p = 806 MeV , (16) 

and A is an ultraviolet cutoff. As one could expect from (|13|) , the cutoff sensitivity of ([TB]) 
is determined by the value of the combination (1 + 2a + (3). Indeed, for large values of the 
loop-momentum, the integrand in (O) is proportional to 



f(f,f) £2 ^l + 2a + p. (17) 

The following subsections are devoted to the determination of a and (3. At first we shall 
analyse the experimental information coming from Kl — > , ~y£ + £~ and Kl — > p + p~e + e~. 
Then we will constraint the value of (1 + 2a + (3) by studying the behaviour of f(q 2 , q 2 ) 
at large q 2 in the framework of perturbative QCD. Finally we will discuss the consistency 
of the previous findings with a model-dependent determination of a and (3 within the 
approach proposed in |2"0f. 
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3.1 Experimental determination of a and j3 



As anticipated, a and (3 are in principle accessible by experiments in the decays Kl —> 
and Kl — > fi + fi~e + e~, dominated by the ifx — >■ 77* and — > 7*7* form fac- 
tors respectively The differential decay rates of K L — > •y£ + £~ and K L — > fi + fi~e + e~, 
normalized to Tj 7 = T(Kl — > 77), are given by 

1 <iri " 2 ^"""^/(^O^A^flX.o) G £ (g 2 ), (18) 
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7e + e~ [|2l], HH and ^ — > 7/i + /i~ |23|] let us extract 
useful information about the q 2 dependence of f(q 2 ,0). The experimental results have 
been analysed up to now assuming only the form factor proposed by Bergstrom, Masso 
and Singer (BMS model) J7|. The latter depends on one unknown parameter a* K and, 
expanding in powers of q 2 /m 2 , can be written as 



f(q 2 , W ^ 1 + (1 " 3.1«^^ + O (tf)*fm\ 

mi v 



(22) 



The fitted values of a* K are given by 



a 



K 



— 0.280 ± 0.083 
— 0.28 ± 0.13 £2 1 
-0.028 l^m IH 



+0.054 
-0.034 



[21] 



(23) 



and the corresponding weighted average is 

a* K = -0.204 ±0.062 . 



with the one proposed in 



(24) 

we obtain the fol- 



(25) 



Comparing the BMS form factor 
lowing relation 

a= -1 + (3.1 ± 0.5)a# , 

where the error is due to the different quadratic dependence on q 2 /m 2 . Then, using (]2 
we find 

a = -1.63 ±0.22 . (26) 
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As already pointed out in |[20| , it must be stressed that an improved determination of 
a would be possible if present data were not analysed assuming only the BMS model. 
Furthermore, an experimental determination of the quadratic term in the expansion of 
f(q 2 , 0) would be extremely useful in order to perform a consistency check of our approach. 

Contrary to a, the experimental determination of f3 is much difficult. In principle the 
K L — > e + e~fi + fi~ rate should be sensitive, in the region where both di-lepton pairs have 
a large invariant mass, to the higher structure in momenta carried by the (3 component 
of the form factor. However, the real sensitivity of this process to (3 is rather small. Thus, 
even though the first evidence for Kl — > e + e~fi + [i~ has been recently reported it is 
unlikely that f3 will be measured with a reasonable accuracy in the short term. 



3.2 Perturbative evaluation of f(q 2 ,q 2 ) 

In the limit q\ = q\ = q 2 3> m 2 K we can simply evaluate the form factor within perturba- 
tive QCD. At the lowest order in a s , the only diagrams that contribute to f(q 2 ,q 2 ) are 
those shown in Fig. 2 ||. Neglecting masses and momenta of the external quarks, as well 
as the contribution of the top quark inside the loop (suppressed by CKM factors [fLlp ), 
the result can be written as 



f QCD (q 2 ,q 2 )=N, 



9a 



Am 2 :. 



Am 2 



(27) 



where 



and 



9g{r) 



r—J(r) + 
dr 



l + 2r 1 M 2 , 

— J{r) + -ln^r 

or 3 m 



2 J 1/r — 1 arctan Jr/(1 — r) + 2 



, 1 - Jl - 1/r 
1 - 1/r In \ +iir | +2 



J(r) 

Jl - 1/r f In 

1 + y/l -1/r 
The normalization factor of (|27]) is given by 



(2* 



< r < 1 



r > 1 



(29) 



16 XGpF^aem 
"9"|F(0,0)|7rv^ 



0.20 



(30) 



where A denotes the sine of the Cabibbo angle JT2[ and F n ~ 93 MeV the pion decay 
constant. The first term in ( p8[) is the contribution of the diagram in Fig. 2a, whereas the 
second one is originated by the graphs in Fig. 2b-c. We have neglected all the contributions 
independent of quark masses that cancel via the GIM mechanism and, whenever possible, 
we have considered the limit My/ — > oo (this is always possible except for the \n(Myy/m 2 ) 
term originated by the diagrams in Fig. 2b-c). 
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Figure 2: Lowest-order quark diagrams that contribute to the Kl — > 7*7* transition 
(every diagram is also understood with the corresponding photons crossed). 



From the above equations it follows 



0{m 2 Jq 2 ) q 2 > 4m 2 

\^ef QCD (q 2 ,q 2 )\ = \N F \l 14 1 , m 2 A „ 2 A . (31) 

1 1 — + - In 4m 2 < g 2 < 4m 2 . 

9 3 <^ 



Using this approximate expression in (|13|) and keeping in the final result only the dominant 
ln(m 2 /m 2 ) terms, leads to the approximate formula of Voloshin and Shabalin for $teAi ong 
II . This result indicates that the long-distance dispersive amplitude of Kl — * is 
very small, however it cannot be trusted in detail since the low q 2 limit of f® CD {q 2 ,q 2 ) 
is completely out of control in perturbative QCD. A more detailed analysis of ^keAi on g 
at the quark level has been recently presented in ||11|| , where the leading QCD correction 
has been estimated. Nonetheless, also the final result of [3 1] cannot be considered fully 
conclusive since an arbitrary infrared cutoff is introduced in order to avoid the dangerous 
low q 2 region. 

As anticipated, our strategy is to use /** (q 2 , q 2 ) to fix the high q 2 behaviour of the 



low-energy parameterization (|T4"D. The simplest requirement that we can derive from (|27| ) 
is that f(q 2 , q 2 ) must vanish for q 2 > 4m 2 . This condition can be implemented in the 
phenomenological expression ( |I"5D in two ways: in a weak sense, assuming 

A 2 < 4m 2 , (32) 

or in a strong one, imposing the "sum-rule" 

1 + 2a + (3 = . (33) 



To be conservative we will use only the weak bound in (p2D , the strong one would have 
been correct only if the low energy parameterization ( fT4] ) was valid also above the charm 
threshold. A more realistic constraint on |1 + 2a + /3\ can be obtained imposing the 



matching between fll4|) and (j27j) for Aq CD <g < 4m 2 . In this case from the second line 
of ( j31~D we obtain 

14 

\l + 2a + f3\ ~ — \N F \ -0.3 . (34) 
9 
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Interestingly, this result suggests that the sum-rule (|33|) is violated only in a mild way 
below the charm threshold. We recall, for comparison, that naive dimensional arguments 
could not exclude values of |1 + 2a + f3\ one order of magnitude larger than in fl3"4"|). The 
smallness of |1 + 2a + /3| is further supported by the leading QCD correction to fQ° D . 



Indeed, as discussed in || [TT[| , the main effect of this correction is an overall multiplicative 
factor smaller than one. 

Combining (|32|) and fl3~4]), we believe that a realistic bound for the last term in (|1^) is 
given by 

|l + 2a + / 3|ln(A/m /3 ) < 0.4 . (35) 

We finally note that is not possible to fix the absolute sign of (l + 2a + (3) in the framework 
of perturbative QCD. Indeed, since we cannot trust the low q 2 limit of the perturbative 
calculation, we are not able to fix the relative sign between the unnormalized form factor 
(F QCD (q 2 , q 2 )) and the K L — > 77 amplitude {F(0, 0)). Moreover, the sign of F(0, 0) is not 
clear also in the framework of Chiral Perturbation Theory due to the cancellation of the 



lowest-order contributions to A(Kl — > 77) [19|, [2 



3.3 Determination of a and j3 in the FMV model 

A more precise, but also more model-dependent, determination of a and (3 can be achieved 
within specific hadronization models. The Factorization Model in the Vector couplings 
(FMV) was proposed in |2(| as a framework to compute the factorizable contributions to 
weak vertices involving vector mesons. This model was proven to be efficient in achieving 
a satisfactory joint description of the vector meson exchange contributions to K — > ^77 
and K L — > 77*, giving 



<*fmv = - jj= °\p^ (2A V + lv)fW - -1-22 (36) 
quite close to the experimental result in fl26|) f\ In (^) £ v = Sfyrriy /(16\/2ii 2 F 2 ), f v is 



fixed from r(p° — > e + e ) to be \ fy\ — 0.20, T(uj — > n 0r y) gives \hy\ — 0.037 and m v 
Moreover, G 8 ~ 9.2 x 10 -6 GeV -2 is the effective coupling of the octet 0{p 2 ) weak chiral 



Lagrangian determined from K — > ixix and rj ~ 0.21 was fixed in p(J from the weak VP^y 
vertex. The application of this model to the construction of the Ki 7*7* vertex through 
vector meson dominance and, consequently, the Pseudoscalar- Vector- Vector (PVV) weak 
vertex is straightforward and gives (assuming only octet contributions) 

256vr G 8 a em m 2 v 3 

t3FMv = jjiFAFm\ SvhvvlAi - (37) 



b In |2(]] it was concluded that the proper estimate of the Kl — » j( + £~ slope should be obtained 
adding to the FMV prediction a contribution generated by a weak Vector- Vector transition (the main 
ingredient of the BMS model) since they are independent contributions. However the two terms have a 
different momentum structure and the BMS one is negligible at large q 2 , i.e. in the region where we are 
interested in the value of (I + 2a + (3). 
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Note that the experimental value of the 7r° — ► 77* slope implies fyhy > 0, thus the relative 
signs of olfmv-, Pfmv an d A(K L — > 77) are completely fixed. The overall arbitrary sign 
is constrained by the experimental data on a (i.e. imposing o>fmv to be negative). 
Combining the predictions of a and (3 in the FMV model we get 

1 + 2a FM v + Pfmv = -0.01 . (38) 
This result is perfectly consistent with the QCD bound in (]35f) . 



4 Numerical results 



The theoretical bound on (1 + 2a + 0) \n(A/mp) in (p5|) , together with the experimental 
determination of a in (|2"6|), let us estimate 1 3?e^4/ on9 1 by means of ( |Tj ) and (|TJ). In order 
to combine the two pieces of information we must assume a statistical distribution for 
(l + 2a + /3) h^A/m^). Assuming for the latter a flat distribution between —0.4 and +0.4, 
and combining it with the Gaussian distribution of a, we find 



long 



< 2.9 x 10" 



(90% c.l.; 



(39) 



The same result is obtained assuming for (1 + 2a + (3) ln(A/mp) a Gaussian distribution 
with central value and a = 0.8/Vl~2 (the a of the original flat distribution). However, 
in this case one can distinguish better the various contributions to the limit fl3"9"|). Indeed 
we find 



\$le A, 



long I 



2al m ml(3,B{K L - 77) 



1 1/2 



2 rn 2 



K 



5.25 + 3.47a + 3(1 + 2a + 0) In 



A 



(40) 



1.61 x 10" 5 x 10.41 ± 0.76 ± 0.691 = 10.66 ± 1.651 x 10" 



Interestingly enough, the knowledge of the absolute sign of the central value of ^R,eAi on g 
(i.e. the relative sign between short and long distance contributions) is not very important 
at this stage, given the large value of the error in fl4"0). Moreover, the above expression 
shows that at present the largest source of uncertainty is generated by the experimental 



error on a. 



Having derived a numerical estimate for |3?e^4; ori9 | we are finally able to extract some 
short distance information from the measured value of B(Kl — > yU + /i~): 



1. Bound on p. Using the Bayesian prescription of the Particle Data Group [ 17| , we 
construct a statistical distribution for |9fte^4 e2: p| that eliminates the unphysical val- 
ues. Then, combining it with the Gaussian distribution of $t,eAi ong discussed above, 
we obtain a distribution function for (|3fJe^4 ea; p| + \$t,eAi on g\) ■ Finally, using this dis- 
tribution in (0), together with m t {m t ) = (167 ±6) GeV and |V c6 | = (0.040 ± 0.003) 



14 , we find 



p > -0.38 



or 



p > -0.42 (90% C.L.) 



(41) 
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This result is consistent with the recent analysis of the CKM matrix presented in 
14"]. Unfortunately the combined constraints on p coming from 6%, ^v(iB d and \V u b\ 



already indicate p > —0.4. However, it must be stressed that the bound in (f|l|) is a 
90% C.L. limit and it has a different statistical significance than the one obtained 
in 



14] scanning on ±<r intervals of the input values. 



2. Bound on \$teUd s \- Similarly to the previous case we can derive a bound on I'SteUds] 
by means of Eq. (^). Treating also l^cA^I as a statistical variable (assuming a flat 
distribution for p between —1 and +1) we find 

\$leU ds \ < 2.4 x 10" 5 (90% C.L.) , (42) 

confirming the original result of Nir and Silverman || . 

To conclude, we stress that the bounds in (O) and ( f£2D must be considered only as a 
preliminary result which can be improved both on the theoretical and experimental sides. 
Indeed, as discussed above, the largest source of uncertainty in the estimate of |9fte^4i „ 9 | 
is given at present by the error on a. Moreover, a bound on (3 and/or a measurement of 
the quadratic slope in K L — > •y£ + £~ could provide interesting consistency checks of our 
approach. Last, but not least, a more stringent bound on |^fte^4 ea; p| could reduce the error 
in the extraction of short-distance parameters. Thus a substantial improvement could be 
foreseen in the near future with the advent of new high-precision data in the sector of 
rare K L decays. 
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